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Using a recently developed strong-coupling method, we present a comprehensive theory for dou- 
blon production processes in modulation spectroscopy of a system of ultracold fermionic atoms in 
an optical lattice. The theoretical predictions compare very well to the experimental time traces 
of doublon production. For experimentally feasible conditions, we provide a quantitative prediction 
for the presence of a nonlinear "two-photon" excitation at strong modulation amplitudes. 

PACS numbers: 03.75.-b, 03.75.Ss, 71.10.Fd 



I 



I 

o 



> 

O 
oo 

^' 

O 
en 

> 



Introduction. Lattice modulation spectroscopy has ad- 
vanced to a standard technique in the physics of ultracold 
atoms in an optical lattice [THT7] . Of particular interest 
is the possibility to measure the Mott gap and the cre- 
ation and analysis of long-lived doublons in a fermionic 
Hubbard model with strong repulsion [^16). 

On the theoretical side, DMRG work has been done in 
ID systems [7. Based on Fermi's golden rule, it has also 
been possible to utilize an equilibrium theory to estimate 
correct doublon production rates [TSl [H]. Nonequilib- 
rium dynamical mean-field theory calculations have ana- 
lyzed models with features analogous to lattice modula- 
tion spectroscopy by including time-dependent hopping 
and time-dependent interactions and showing how they 
affect the double occupancy [8j. In the linear-response 
limit, quantum Monte-Carlo calculations have studied 
phase correlations between the double occupancy and the 
lattice modulation [12]. 

In a Mott insulator, doubly occupied sites can be inter- 
preted as occupied by "doublon" quasi-particles. Their 
long life time is due to a separation of energy scales which 
requires a rather rare high-order (in the hopping) many- 
body process for a decay to occur IIJ. In this Letter, 
we provide a computational study of the creation of dou- 
blons due to lattice depth modulation. We derive time- 
dependent tight-binding parameters for a modulated lat- 
tice and then apply a recently developed computational 
strong-coupling method (TB]. We then validate our ap- 
proach by making contact with experimental data by 
Greif et al. [13] who have provided a detailed measure- 
ment of the time evolution of the creation process of dou- 
blons in a ^°K system. Exploring the parameter space 
further, we find that for sufficiently high modulation am- 



plitudes, processes involving the nonlinear combination 
of two coherent quanta of the many-body system en- 
hanced doublon production rates at a frequency which 
equals precisely half the value of the Hubbard repulsion. 
Higher order nonlinear effects are difficult to produce due 
to the way the amplitude modulation of the potential 
translates into the time dependence of the microscopic 
parameters of the single-band Hubbard model. 

Method. The time dependence of the lattice depth 
modulation is set as follows: 



V^W - K) + X[o,t„„a] W ■ AF sinc^i. 



(1) 



where Vo, AV, and u are the average value, modulation 
amplitude, and modulation frequency of the optical lat- 
tice potential depth, respectively, and 



X[o,t„od](0 



1, iffe[0,tn 

0, otherwise 



(2) 



is the characteristic function of the time interval over 
which the lattice depth is modulated. The modulation 
time period length imod = "-mod • Stt/oj is a function of the 
number of modulation cycles nmod chosen for the driving 
of the system. 

The Hamiltonian for a single atom in a d-dimensional 
optical lattice is given by [TS] 

i^.ingicW - -7^^^ + V{t) • Vsin2(A:x,), (3) 

where k = 27r/A is the lattice vector, with the laser wave- 
length A. A natural energy unit to use is the recoil-energy 
Efi = fi^k'^/2m. Using the respective time-dependent 



maximally localized Wannier functions [20], we map the 
Hamiltonian to a single-band lattice model. Note that for 
certain frequencies and amplitudes, transitions to higher 
bands will eventually become important. For inter-band 
transitions in particular, corrections from terms coming 
from the time-derivative of the Wannier functions have 
to be taken into account [T7]. But this is not needed for 
the case we evaluate here, as we always keep the system 
in the single-band limit. 

The many-body physics of fermionic atoms with spin 
1/2 is then described by the single-band fermionic Hub- 
bard model [H] 
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H{t)=-J{t) J2 (clc,.+h-c. 



U{t) ^ ni^n^i + ^ e{t)7ii„, 



(4) 



where J{t) = -{wi{t)\Hsingic{t)\wj{t)) is the hopping 
between Wannier states Wi{t) and Wj{t) at neighbor- 
ing sites i and j, e(t) — e(t) — U{t)/2 — fi is the 
on-site energy with e(i) — {w i{t)\H sing\c{t)\w i{t)) , and 
U{t) = g j A''-r\wi{r;t)\'^ is the time-dependent repulsion 
of atoms, while g — Airh^a/m is determined by the s- 
wave scattering length a ^. The bracket (•,•) denotes 
nearest-neighbor pairs on the lattice. Note that the as- 
sumed time-dependence of the single-particle energies is 
unimportant, since in this work, we perform calculations 
on a homogeneous, translationally invariant lattice with 
a fixed density of fermions, and in this case equal-time 
expectation values are independent of the time evolution 
of e(i). The comparison to experiment is made within 
the Mott insulating region, where we set the chemical 
potential to /j, = e(0), for half filling. 

In order to compute non-equilibrium observables as a 
function of time we use a strong-coupling approach which 
self-consistently expands the self-energy to second order 
in the hopping .18 j . The formalism enables us to numer- 
ically evaluate the on-site contour-ordered Green's func- 
tion 



GAt,t') = -i{Tcc,{t)cl{t')), 



(5) 



where t and t' are times located on the Kadanoff-Baym- 
Keldysh contour C 18J. In order to evaluate the site's 
double occupancy D{t) — {n^n^){t), we use the following 
relation for its equal-time derivative, where t <c t': 



dGAt,t') 



dt 



U{t)D{t)+e{t){n,){t) 



^kin 



{t). (6) 



t'=t+ 



In this expression, the contribution of the spin state a to 
the kinetic energy per atom can be evaluated in momen- 
tum space via 



s^-(i) = -2-^^(nfe.)(t) Vcosfc„ 

'^ k 



' /^' 



(7) 
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FIG. 1: (color online) Comparison of the double occupancy 
as a function of time to the experimental data from Ref. [13) . 
Given parameters of the experiment are U0/6J0 = 4.1, Vb = 
7Ej{, Ay/Vo = 20%, hw = Uo- Since the initial temperature 
T is the only unknown parameter, we plotted results of our 
theory for several values of T. The temperature decreases 
from top to bottom as indicated by the labels. Error bars 
denote statistical errors from multiple measurements. 



Comparison to Experiment. In order to validate the 
strong-coupling approach. Fig. [l] provides a comparison 
to recent experimental data on the process of doublon 
creation due to lattice modulation in a 3D Hubbard 
model with '^''K [T3]. The figure shows the double occu- 
pancy D{t) as a function of time as the lattice depth V{t) 
is modulated. The initial value of V{t) is Vq = 7-Efl, and 
it is modulated by Al^/Vb = 20% at a frequency huj = Uq. 
The two-body scattering length is tuned through a Fes- 
hbach resonance such that U0/6J0 ~ 4.1. For an ini- 
tial temperature fc^T = 0.15Uq, the double occupancy 
matches the initial experimental value at i = 0, and 
the theory yields a good description of the subsequent 
time-dependence. Since the lattice modulation frequency 
equals the Hubbard repulsion Uq — U{0), particles are 
resonantly excited from the lower to upper Hubbard band 
during the lattice modulation. There also exists a process 
of de-excitation which is represented by the decreasing 
sections of the observed curve. For the system studied in 
Fig. [1] the excitation process dominates. 

Fig. [2] shows the same analysis for different modula- 
tion frequencies. Clearly, the doublon production rates 
at these frequencies are much lower than for the resonant 
case hiu — Uq. As a consequence, the signal-to-noise ra- 
tio in the experiment is also higher. We again find that 
fcsT w 0.15t/o is a reasonable estimate for the initial 
temperature. In some cases, however, the data seem to 
be more compatible with fc^T = O.IUq or ksT = O.OTC/q. 
Although the experimental preparation procedure of the 
initial thermal state was identical for each data point. 
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FIG. 2: (color online) Same as Fig. [I] for off-resonant mod- 
ulation frequencies. Data provided by the Esslinger group. 
The temperature decreases from top to bottom in each panel. 
Error bars here are just the spreads from three successive 
measurements. 
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FIG. 3: (color online) Effect of the amplitude on lattice pa- 
rameters and doublon production. Model parameters are 
Vb = 15-Efl, C/0/6J0 = 31.10, ksT = 0.2(7o, hw = [/q. The 
modulation amplitude is either 25% or 50%, as specified for 
each graph. 



systematic drifts in the temperature are quite possible 
and hard to predict or measure. 

Amplitude Effects. Let us now discuss the effect of 
the amplitude ISV in more detail. In order to pre- 
vent inter-band transitions from becoming important, 
we study this problem for a deeper lattice, such that 
mintg[o,t^,^j] V{t) > 7Eji. We have plotted the lattice 
depth V{t), the normalized hopping J{t)/U{t), and the 
double occupancy D(t) as a function of time for a deeper 
lattice at two strong values of the lattice modulation am- 



FIG. 4: (color online) Double occupancy as a function of time 
for lattice modulations at different frequencies. The other 
model parameters are same as in Fig. 3] for AV^/Vb = 50%. 
Vertical bars denote the respective end time of the modula- 
tion. The modulation frequency associated with a given graph 
is shown with a label on the curve. 



plitude in Fig. [3j As the amplitude is raised, the non- 
linear relationship between the hopping and the lattice 
depth results in a periodically kicked rather than a peri- 
odically driven system (due to J(t)/U{t) becoming very 
small for deep lattices). Overall, the double occupancy is 
increased stepwise within each modulation cycle. Again, 
we observe both, excitation and de-excitation processes 
in the double occupancy data. Within the first mod- 
ulation cycle, we observe the former and the latter to 
coincide with the decreasing and increasing regimes of 
J{t)/U{t), respectively. Within modulation cycles fur- 
ther out in time, the stepwise increase in double occu- 
pancy starts already when J[t)/U{t) is still rising. This 
is due to the fact that between the spikes in J{t)/U{t), 
the hopping is effectively zero, so that the system oscil- 
lates internally with frequency U{t) sa f/g — f^ within 
the four-level Hilbert space associated with a single lat- 
tice site [32]. As a consequence, there is a constructive 
interference of the immediate effect of lattice modulation 
and the internal oscillation, once the latter is fully en- 
gaged. 

Amplitude versus Frequency. We would like to address 
the interplay of the internal oscillation at zero hopping 
and the lattice modulation in more detail next. Let us 
discuss data at different frequencies by first assuming a 
strong lattice modulation amplitude AV/Vq — 50%. The 
modulation time interval in Eq. ( ^ for a given frequency 



is chosen to be t 



mod 



• — . Results are shown 
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in Fig. |4] The strongest increase in double occupancy is 
observed in the resonant case Huj = Uq. As the end time 
^mod of the modulation is exceeded, the system contin- 
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FIG. 5: (color online) Doublon production as a function of 
frequency for different amplitudes. Amplitude increases from 
bottom to top. The other parameters are the same as in 
Figs. [3] and m 



modulation spectroscopy due to the non-linearity of the 
map V — >■ J/U. Its presence could be possibly enhanced 
by designing V{t) in such a way that a harmonic shape 
is obtained for J{t). 

Summary. We have studied finite-amplitude lattice 
depth modulation spectroscopy of ultracold fermionic 
atoms in the Mott-insulating phase using a recently de- 
veloped strong-coupling method [TB] ■ In order to validate 
the theory, we have compared to experimental data and 
found excellent agreement. Only the temperature of the 
initial thermal state was unknown in the experiment and 
had to be determined a posteriori. We furthermore an- 
alyzed higher amplitudes of the modulation strength for 
deeper lattices. A large value of the amplitude results 
in a pulsed, rather than a driven system, in terms of 
the time-dependence of the hopping. This causes step- 
shaped changes in the double occupancy, accompanied 
only by oscillations with frequency U{t) ~ Uq of local 
degrees of freedom on a single lattice site. At a certain 
threshold in amplitude, a second " nonlinear" peak in the 
doublon production rate appears at cj = Uo/2. 



ues oscillating with the internal frequency Uq. A signifi- 
cant increase in double occupancy is also observed in the 
case hoj = Uo/2. Moreover, in this case also the sys- 
tem continues oscillating with C/g after modulation, but 
with a smaller amplitude. For the off-resonant frequen- 
cies Huj — 1.43Uo and Huj — 0.77?7o, some intermediate 
excitations to the upper Hubbard band show up but (par- 
tially) annihilate after a while. The amplitude of internal 
oscillations after the modulation is turned off appears to 
grow with the total increase in double occupancy reached 
when the modulation ends. 

To provide a further overview of the frequency and 
amplitude dependence. Fig. [5] shows the final value of 
the double occupancy, i.e. the value averaged over one 
oscillation cycle [imodjimod + h/Uo], vs. the frequency 
for different amplitudes as a lattice modulation "spec- 
troscopy" . Evidently, a second-order resonance is ob- 
served at frequency huj = O.bUo, which is suppressed as 
the amplitude is lowered. It can be interpreted as a coher- 
ent excitation involving the nonlinear combination of two 
smaller quanta of energy Uo/2, i.e. an analogue to two- 
photon excitations in quantum optics. As a consequence 
of the combination of these two quanta, the width of 
the second-order peak is approximately half of the width 
of the first-order peak. Also note that both peaks are 
shifted towards slightly smaller frequencies due to the 
finite width of the Hubbard bands. Lower-energy excita- 
tions are possible by exciting from the upper edge of the 
lower Hubbard band to the lower edge of upper Hubbard 
band. 

Whereas the Uo/2 peak is clearly visible for very large 
modulation amplitudes, the rather isolated time evolu- 
tion of lattice sites between the kicks in J{t)/U(t) pre- 
sumably suppresses its amplitude. This is inherent to 
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